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Abstract 

Here homogenization theory is used to establish a connection between the symmetries 
of a periodic elastic structure associated with the microscopic properties of an elastic 
material and the material symmetries of the effective, macroscopic elasticity tensor. 
Previous results of this type exist but here more general symmetries on the microscale 
are considered. Using an explicit example, we show that it is possible for a material to 
be fully anisotropic on the microscale and yet the symmetry group on the macroscale 
can contain elements other than plus or minus the identity. Another example demon¬ 
strates that not all material symmetries of the macroscopic elastic tensor are generated 
by symmetries of the periodic elastic structure. 

Keywords: microstructure, macroscopic elasticity tensor, material symmetry, multi¬ 
scale analysis 

MSC subject classification: 74B05, 74Q15 


1 Introduction 

Elastic composites span a wide range of materials, from geometerials and biological tissues 
to synthetic materials. See, for example, Cherkaev and Kohn [3] and Jones HU. Knowing 
the properties of such materials can lead to a better understanding of their behavior in 
their physical environment. The properties of a composite are dictated by the properties 
of its constituents. However, determining the properties of a composite from those of 
its constituents is a nontrivial matter. Homogenization theory can be used to rigorously 
derive the macroscopic properties of a composite from the microscopic properties of its 
building blocks. Here we are interested in the case where the composite consists of a 
periodic arrangement of linearly elastic materials. In this instance, homogenization theory 
provides techniques for taking the limit as the size of the periodic microstructure goes to zero 
and obtain an macroscopic elasticity tensor associated with the composite. Determining 
the macroscopic or effective elasticity tensor involves solving systems of elliptic partial 
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differential equations on a representative unit cell domain associated with the periodic 
microstructure of the composite material. See, for example, Cioranescu and Donato |1], 
Jikov, Kozlov, and Oleinik |10) . and Oleinik, Shamaev, and Yosifian m for results on the 
homogenization of the equations of linear elasticity. 

Besides having a microscopic structure, some important engineering and biological mate¬ 
rials have symmetry properties. It is possible for the response of a material (not necessarily 
a composite) to be unaffected by a change in reference configuration. This leads to the 
notion of material symmetry. Roughly speaking, the material symmetry group consists of 
all transformations of the reference configuration that leave the response of the material 
unchanged. A connection between material symmetry and microstructure is well-known as 
different material symmetry groups are connected with different crystalline structures. For 
a list of which material symmetry group is appropriate to assume for a given crystalline 
structure see, for example, Coleman and Noll [5] or Curtin j^. Thus, it is reasonable 
to conjecture, in the context of homogenization theory, that the properties of a mate¬ 
rial on the microscopic scale yield information about the material symmetry group on the 
macroscale. In this paper we show that if the elasticity tensor describing the properties 
of the microstructure satisfies an invariance condition involving an affine transformation 
that preserves volume, see equation (|^ in Section then the gradient of this transforma¬ 
tion is a material symmetry of the macroscopic elasticity tensor. Previous results of this 
type were established by Jikov, Kozlov, and Oleinik |in| and Alexanderian, Rathinam, and 
Rostamian [1], but these authors only considered transformations that are rotations or re¬ 
flections about a fixed point. Examples of transformations not considered before include a 
translation together with a rotation or a unimodular transformation that is not a rotation. 
A detailed discussion of how previous results compare with what is established in this paper 
is given after Theorem |3.2| in Section 

One of the advantages of knowing the material symmetry group of the macroscopic 
elasticity tensor is related to numerics. The components of the macroscopic tensor relative 
to a basis are obtained by solving unit cell problems, which are systems of elliptic partial 
differential equations. Taking into account the major and minor symmetries of the elas¬ 
ticity tensor, the 21 elasticities are determined by solving six unit cell problems. Hence, 
information that reduces the number of unknown components can decrease the number of 
unit cell problems needed to be solved and, thus, save considerable computational time. 
Knowledge of the material symmetry group yields this kind of information. For example, if 
the macroscopic elasticity tensor possesses cubic symmetry, the nine unknown elasticities 
can be computed by solving only two unit cell problems. 

The outline of the paper is as follows. In Section the concept of a periodic elas¬ 
tic structure and an appropriate notion of symmetry are introduced. Section contains 
the main result connecting symmetries of the periodic elastic structure with the material 
symmetry group of the macroscopic elasticity tensor. In Section several examples of 
periodic elastic structures are given and the corresponding material symmetry groups are 
mentioned. One example showns that it is possible for the constituents of a composite to be 
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completely anisotropic and yet the macroscopic elasticity tensor can have a symmetry gronp 
containing elements other than pins or minns the identity. Another example demonstrates 
that there can be elements of the material symmetry gronp of the macroscopic elasticity 
tensor that are not associated with any symmetry of the periodic elastic strnctnre. Finally, 
Section contains some conclnding remarks. For completeness, we inclnde a derivation of 
the macroscopic elasticity tensor in the Appendix. 

2 Periodic elastic structures 

In this section we introduce the concept of a periodic (linearly) elastic strnctnre and discnss 
the appropriate notion of symmetry associated with this strnctnre. Snch symmetries are 
generalizations of the transformations nsed when discnssing material symmetry. 

For the pnrposes of this paper it is nsefnl to distingnish between points and vectors 
and, hence, work in a Enclidean point space £ with a vector space V eqnipped with an 
inner-prodnct. Lowercase letters snch as y or ^ will denote elements of £ and are points, 
while boldface lowercase letters snch as a and b denote vectors. Boldface nppercase letters, 
snch as H, will denote linear mappings from V to itself. The inner-prodnct on V will be 
nsed to identify V with its dnal. It follows that given a linear mapping H from V to itself, 
we can consider its transpose also to be a linear mapping from V to itself. 

Consider a parallelepiped Y \n £. The cases when £ is two- or three-dimensional are 
of primary interest, however here it is assnmed that £ is n-dimensional, with n G N. 
The set Y can be nsed to form a periodic tessellation of £. Let bj, i = 1, ... ,re, be the 
linearly independent vectors in V that are the length of and parallel to the edges of Y ; see 
Fignrej^a). For all 2 : G T, there are nniqne integers fcj, i = 1, ..., n, snch that 2 — fcibj 
is an element of Y. Thns, we can define 

n 

{zjy := 2 - ^/cjbj G y. (1) 

i=l 

Here we are interested in linearly elastic materials that have a periodic structure. To 
this end, consider a position dependent elasticity tensor C defined on Y and extended 
T-periodically to all of £ so that 

C{z) = C({ 2 ;}y) for all z ^ £. (2) 

The pair (Y, C) is referred to as a periodic (linearly) elastic structure. In such a construction, 
Y is called the unit cell. 

Consider another parallelepiped Y in £ and elasticity tensor C defined on Y so that 
(y, C) is another periodic elastic structure. We say that the periodic elastic structures 
(y, C) and (y, C) are equivalent if the periodic extension of C relative to Y is equal to the 
periodic extension of C relative to Y. With the aid of Q, this can be expressed as 

C({ 2 }y) = C({ 2 }y) for all z G T. 
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It is possible for the two periodic elastic structures to be equivalent even when Y ^Y. To 
see this, consider a periodic elastic structure (Y, C) and let a be a vector. Extend C as in 
([^ and dehne Ca to be its restriction to the set a + T, which is the translation of the set 
Y by a. The structure (a + Y, Ca) is equivalent to (Y, C). 

Here we consider symmetries of a periodic elastic structure generated by affine trans¬ 
formations that preserve volume. Such a transformation h of £ has the form 

h{z) = Zo + sl+ H( 2 ; — Zo) for all z £ £, (3) 

where Zo is a point in T, a is a vector, and H is a unimodular linear mapping, so that 
|detH| = 1. To interpret h, hrst consider the case where H is a rotation. In that case, 
in ([^ a acts as a translation and Zo is the point about which the rotation H takes place. 
When H is not a rotation, a similar interpretation applies but instead of the space being 
rotated about Zo, it is transformed about 2;o by H. Unless otherwise stated, when speaking 
of a symmetry h, it is assumed to have the form 

Some explanation to why transformations of the form Q is warranted. When discussing 
material symmetry of a material at a point Zo, transformations of the form 


g{z) = Zo + G(z - Zo) (4) 

are considered. See, for example, equation (50.6) on page 285 of Gurtin, Fried, and Anand 
[ 3 . It is assumed that G is unimodular so that g preserves volume. This ensures that the 
mass density of the material at Zo is preserved. See page 11 of Epstein and Elzanowski [3- 
The transformations considered in (|^ also include a translation by a vector a since we are 
interested in symmetries of periodic structures, which, by dehnition, are invariant under 
certain translations. One can think of the transformations of the form ([^ as a merger of 
the transformations considered for material symmetry and those for periodic structures. 
Let Sym denote the set of all symmetric linear mappings from V to itself. 

Definition 2.1. A transformation h of the form (|^ is a symmetry of the periodic elastic 
structure (T, C) if 

C(z)E = H[C(h-^(z))(H^EH)]H^ for all z £8, E e Sym. (5) 

To understand the motivation behind Dehnition |2.1[ consider in ([^ instead of h the trans¬ 
formation g dehned in (Q instead of h in ([^. Then (|^ at z = Zo would read 

C(zo)E = G[C(zo)(G^EG)]G^ for all E e Sym, 


which would imply that the material point Zo has G as a material symmetry. See, for 
example, equation (52.21) on page 300 of Gurtin, Fried, and Anand |9]. Thus, Dehnition 


2.1 should be viewed as a generalization of the dehnition of material symmetry to periodic 


elastic structures. 
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(a) 


(b) 

Figure 1: (a) A two-dimensional example of a unit cell Y and the periodic tessellation 
associated with it. The vectors bi and b 2 are the length of and parallel to the sides of Y 
and generate the tessellation in that translating Y by all integer linear combinations of bi 
and b 2 covers the plane, (b) An example of a periodic elastic structure made up of two 
isotropic elastic materials that can be generated by two different unit cells Y and Y that 
are related through h, which is a symmetry of the periodic elastic structure. 


It is easily seen that the collection of symmetries of a periodic elastic structure form a 
group under function composition. For any unimodular linear mapping H, define 

Ch( 2 )E := H[C(z)(H^EH)]H^ for all z G T, E e Sym. (6) 

If C is periodic, then so is Ch- Using the notation in ([^, h is a symmetry of the periodic 
elastic structure if and only if 

C{z) = (z)) for all z ^ £. (7) 

Notice that translations by integer linear combination of bj, where i = 1,... ,n, are sym¬ 
metries of the periodic elastic structure. 

Proposition 2.2. If h is a symmetry of the periodic elastic structure (Y,C), then C is 
hi-periodic where bj := Hbj. 

Proof. Consider z ^ £ and set z := h~^(z). From ([^ and the fact that C is bj-periodic, 
we have 


C(z + b,) = C{h{z + bj)) = Cn{z + bj) = Ch(^) = C(z), 
and so C is bj-periodic. □ 
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The next result says that periodic elastic structures related by symmetries are equiva¬ 
lent. 

Proposition 2.3. Let (T, C) be a periodic elastic structure with symmetry h, and set Y := 
h{Y). If C is defined by 

C{y)-.= Cn{h-\y)) for all y G T, (8) 

then the periodic elastic structures (T, C) and (Y, C) are equivalent. 

Proof. It must be shown that the periodic extensions of C and C relative to Y and Y, 
respectively, are equal. Let z G P be given and set bj := Hb,, for i = 1,..., n. Define {5}y 
analogous to Q, so that there are integers ki, i = 1,... , n, such that 

n 

{z]y = z-'^kihi gY . ( 9 ) 

i=l 


By (|^ and the periodicity of C and C, see ([^, and the definition of Ch in Q , we have 

(C(|) = C({z}y) = C^{h-\{z}y)) = Cn{{h-\{z}y)}Y). (lO) 


Applying h ^ to @ and using the definition of bj yields 


h ^({^}y)=^ ^{z) - '^kihi, 

i=l 


from whence it follows that 


{h-\{z}y)}Y = {h-Hz)}Y- 

Using this last equality in ( |10[ ), the periodicity of Ch, and equation Q results in 

Ciz) = Cu{{h-\z)}Y) = Cnih-\z)) = C{z). 

Since the above equation holds for any z G P, it follows that C = C and, hence, the elastic 
structures are equivalent. □ 


A consequence of Proposition 2.3 is that the elasticity tensors C and C are equal in T, and 


so (Y, C) and (Y, C) are equivalent elastic structures. This fact will be used in the next 
section. 

To illustrate the results of Propositions |2 .2 1 and |2 . 3| consider the tessellation of the plane 
depicted in Figure [^b). Define C so that it equals one isotropic tensor on the white triangles 
and another isotropic tensor on the grey triangles. This elastic structure can be generated 
by the unit cells Y and Y shown in the figure. The transformation h that takes T to T is 
a symmetry of the periodic elastic structure. Moreover, the vectors bi and b 2 associated 
with the T-periodicity of the structure are related to the vectors bi and b 2 associated with 
the T-periodicity of the structure through the gradient of h —that is, bj = Hbj, i = 1,2. 
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3 The macroscopic elasticity tensor 


Here we recall the macroscopic elasticity tensor of a composite material obtained by ho¬ 
mogenizing the equations of linear elasticity associated with a composite with periodic 
microstructure. After this, we prove our main result connecting the symmetries of the 
microstructure with the material symmetries of the macroscopic elasticity tensor. 

Homogenization of the equations of linear elasticity is classical. See, for example, muni 
HU. For the sake of completeness, a sketch of the derivation of the macroscopic equations 
of linear elasticity in a domain with periodic microstructure is given in the Appendix. 

Here we need to recall the formula for the macroscopic elasticity tensor. The starting 
point is to consider a linearly elastic material with a periodic microstructure. We consider 
the periodic elastic microstructure given by 

C^{z) = - - -|- qj for all z ^ £, (11) 

where q is an arbitrary, but fixed, point, e is a small parameter associated with the length 
scale of the microstructure, and {Y, C) is a periodic elastic structure. If W is the result of 
scaling Y hy £ about the point q, then (Y^,C^) is a periodic elastic structure, which can 
be viewed as a microstructure since £ is small; see Figure 

Homogenization theory provides techniques to derive an effective elasticity tensor 
that describes a material whose behavior approximates that of a material with periodic 
microstructure [Y^, C*") for small e. The effective elasticity tensor is given by 

= j- C{y) [E -|- Vy'w^(y)] dy for all E G Sym, (12) 

where the symbol J denotes the average integral anc0 G Wper(E, V) is the unique 
solution of 


/y 


C(y)[E. 


VyW 


\y)] ■ ^y^{y)dy = 0 


for all V G Hp^^{Y, V). 


(13) 


See (34) and 


in the Appendix for more details. Since C is an elasticity tensor, it 


possesses major and minor symmetries so that, in particular, for any linear mapping A, we 
have CA = C[^(A -|- A^)]. Thus, only the symmetric part of has an influence in 

([^ and 

The macroscopic elasticity tensor C® can be given in terms of an equivalent periodic 
elastic structure with a different unit cell. In partic ular, if /i is a symmetry of (T, C) and 
Y := h{Y), then the remark made after Proposition 


2.3 


says that (Y, C) is an equivalent 
elastic structure. The homogenization procedure can be carried out using Y instead of Y 
to find that is also given by 


C°E = j- C{y) [E -|- Vj;w^(y)] dy for all E G Sym, (14) 

^See the Appendix for the definitions of the function spaces V) and Wper{Y,V). 
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Figure 2: A depiction of a periodic elastic structure together with a scaling of it. The unit 
cell V generates the initial structure, and the grey and white regions represent different 
elastic materials. The smaller, transparent structure generated by the unit cell is the 
result of scaling the original structure about the point q by e = 1/2. 


where G Wper(T, V) is the unique solution of 


/y 


C(y)[E + Vj 


w 


V^v(y) dy = 0 for all v G Ffper(h", V) 


(15) 


The following result relates the solutions of the unit cell problems (13) and (151 


Lemma 3.1. The solutions of the unit cell problems (13) and (15) satisfy 


Vy(w“ )(/i (y)) = H V^w^(y)H for almost every y €Y = h{Y). 

Proof. Given v G >Vper(T, V), dehne 

v(y) := H^v(/i(y)) for almost every y &Y. 


(16) 


(17) 


Notice that since v G Wper(T, V), then v G Wper(T, V). Hence, considering (13) with E 
replaced by H^EH results in 

^ C(y) [H^EH + V,(wH^™)(y)] • V,v(y) dy = 0. 

Using the change of variables y = h{y) and the dehnition of v in ( |17[ ) yields 

^C(h-Hy))[H’^EH + Vj,(wH^™)(h-i(y))] ■U^Vy^mdy = 0. (18) 

Since /i is a symmetry of the elastic structure we have that C satishes ([^, and so (18) 
becomes 

^ C(y) [E + H-Tv,(wH^™)(/i-Hy))H-^] • Vj^v(y) dy = 0. 


(19) 





















Notice that H ' Vy(w^' ^(y))H ^ = V^(H ™ o/i ^){y) for almost every 

y €Y, and so (19) can be written as 

^ C(y) [E + . v^v(y) dy = 0. (20) 


Since (20) holds for all v G W’per(E, V) and the solution of (15) is unique, this proves the 
lemma. □ 


Using Lemma 3.1 we can now prove the main result of the paper. 
Theorem 3.2. If h is a symmetry of the periodic elastic structure, then 
C°E = H[C°(H^EH)]H^ for all E G Sym. 


( 21 ) 


Proof. Let E G Sym be given and set Y := h(Y). By (12), the change of variables y = h{y), 
Lemma 3.1 the dehnition of a symmetry of a periodic elastic structure ([^, and ( |14[ ), we 
have 


c 


o^H^EH) =£c{y) [H^EH + V,(wH^™)(y)] dy 

= j^C{h-\y)) [H^EH + Vy{^^^^^){h-\y))\dy 

= I C{h-\m [H^EH + H’^VjiwE(y)H] dy 

i^C(y)[E + V^w^(y)]dyH-T 


'Y 

= H 

= H-HC°E)H- 


□ 


Equation (21) is the condition for H to be in the material symmetry group of Thus, 


Theorem |3.2| says that the gradient of every symmetry of the periodic elastic structure is a 
material symmetry of the macroscopic elasticity tensor. For a detailed discussion of material 
symmetry see, for example, Gurtin, Fried, and Anand or Epstein and Elzanowski |7]. 


Notice that Theorem 3.2 does not say that every element of the material symmetry 
group of C® is generated by a symmetry of the periodic elastic structure. In fact, there are 
material symmetries that are not generated by symmetries of the periodic elastic structure. 
This will be demonstrated with an example in the next section. 

Although the above analysis was carried out in the context of linear elasticity, similar 
calculations hold for any hrst-order transport law. In that case, a mobility tensor M would 
be dehned on a unit cell and ([^ would be replaced by 


M(z) = HM(/i-^(z))H 


rT 


for all z G <5. 
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Figure 3: Examples of periodic elastic structures containing symmetries that are not covered 
by previous results similar to Theorem 3.2 In these examples, the white and grey regions 
should be viewed as different isotropic materials, (a) An example of a periodic elastic 
structure having a symmetry consisting of the translation a together with a reflection in 
which neither the translation nor the reflection by themselves symmetries of the elastic 
structure, (b) This periodic elastic structure has amongst its symmetries reflections about 
the vertical and horizontal lines passing through p and rotations about the point q by 7r/2 
radians. 


For example, the physical processes of diffusion, heat conduction, and dielectric induction 
have this kind of structure. 


Results similar to Theorem 3.2 were established by Jikov, Kozlov, and Oleinik m and 
by Alexanderian, Rathinam, and Rostamian [T]. To accurately compare these results with 
take the Euclidean space £ to be M”, in which case points can be considered 


Theorem 3.2 


as vectors and vice versa, 
of the form 


Jikov, Kozlov, and Oleinik proved Theorem 3.2 for symmetries 
z I— )■ Qz for all z G M^, (22) 


where Q is an orthogonal linear mapping and Alexanderian, Rathinam, and Rostamian also 
considered symmetries of this form in the context of the homogenization of diffusive random 
media. Thus, these previous results only considered symmetries that are rotations or re¬ 
flections about a single point—in particular, the origin. Whereas Theorem 3.2 also includes 


symmetries about different points, symmetries consisting of an orthogonal linear mappping 
together with a translation, and symmetries involving unimodular linear mappings. 

A periodic elastic structure with a symmetry consisting of a translation and reflection 
is depicted in Figure j^a). In this figure consider the white and grey regions to consist of 
different isotropic elastic materials. Notice that a reflection about any vertical line is not a 
symmetry of the elastic structure, however the transformation consisting of the translation 
a together with the reflection about the vertical line aligned with a is a symmetry of the 
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periodic elastic structure. Figure [^b) is an example of a periodic elastic structure that 
contains symmetries of the form 

z >->■ Zo + H(z — Zo) for all z S £ 

for different H associated with different points Zo- Namely, there are reflection symmetries 
about point p and a rotation about the point q by 7 r /2 radians is a symmetry. An example 
of a periodic elastic structure that has a symmetry where H is not orthogonal is described 
in the next section. 

When C takes only a hnite number of values, the elastic structure is made up of a 
hnite number of materially uniform constituents. In this case, two factors determine the 
symmetries of such an elastic structure and, hence, the material symmetry of C®: the 
material symmetries of the constituents and their arrangement within the periodic structure. 
This is illustrated in the next section. Parnell and Abrahams |15) realized this in the context 
of a three dimensional elastic structure made of two constituents. 

4 Examples 

In this section several examples of periodic elastic structures are given. For each example, 
the symmetries of the periodic elastic structure are mentioned together with the resulting 
material symmetries of the macroscopic elasticity tensor guaranteed by Theorem |3.2| 

For simplicity, consider T to be the three-dimensional space so that points and 
vectors are interchangeable and use the standard basis 01,02,63. Given a unit vector a 
and an angle (j), let Rf be the right-handed rotation through the angle (j) about the axis in 
the direction of a. Notice that —is the reflection with respect to the plane through the 
origin with normal a. 

Since the macroscopic elasticity tensor is given by an integral expression, if C is 
changed on a set of measure zero, then the macroscopic elasticity tensor remains unchanged. 
For this reason, in the following examples the microscopic elasticity tensor C is only speci- 
hed up to a set of measure zero. Similarly, the unit cell Y which dehnes the periodic elastic 
structure only tiles up to a set of measure zero. 

Orthotropic symmetry: It is possible for the elasticity tensor C associated with the 
periodic elastic structure to be fully anisotropic at each point and yet the macroscopic 
elasticity tensor can have orthotropic symmetry, meaning that 

C°(E) = H[C°(HTeH)]hT for all E G Sym and H = -R;: 3 . 

To obtain this, the anisotropies of the elastic structure are arranged in a particular way. 
Begin by considering an anisotropic elasticity tensor and the unit cell Y = (—1,1)^. 
To dehne C on Y, start by setting C equal to C“ on the positive octant (0,1)^. Then, 
recalling (|^, set C equal to on (— 1 ,0) x (0,1) x (0,1). Next, dehne C{yi,y 2 , 2 / 3 ) with 
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Figure 4: A depiction of a cross-section of the unit cell described in the example involving 
orthotropic symmetry. The cross-section is taken parallel to one of the coordinate axes. The 
anisotropy of the elasticity tensor, depicted as crisscrossed lines, in the different quadrants 
are related through reflections. 

yi G (-1,0) U (0,1), y2 G (-1,0), and 2/3 G (0,1) by C(yi, 2 / 2 , 2 / 3 ) = C_rj ( 2 / 1 ,- 222 , l/s)- 
Finally, define C( 2 /i, 2 / 2 , ys) for negative 223 -coordinate by requiring that C( 22 i, 222 , ys) = 
C-R’T (yi, y 2 , — 223 )- Since reflections about perpendicular planes commute, it is readily 
shown that the mappings 

2 ; !-)■ —for all z G M^, i = 1, 2, 3 

are symmetries of the periodic elastic structure (T, C). See Figure]^ for a depiction of the 
microstructure associated with C. It follows from Theorem |3.2| that —Rg., for i = 1,2,3, 
are material symmetries of C^. 

Tetragonal symmetry: Suppose the unit cell Y = (—1,1)^ consists of two different 
isotropic materials, one representing fibers and the other a matrix in which the fibers are 
embedded. Suppose that there is one fiber that occupies a cylinder whose axis is the y 3 -axis; 
see Figure l^a). Notice that the mappings 

for all z G H = R^^ , Re^^ 

are symmetries of the periodic elastic structure. To have R^^ and Reg^ as material sym¬ 
metries is sufficient for tetragonal symmetry. 

The same symmetry persists if there are several fibers appearing in the unit cell arranged 
as depicted in Figure |^b). In general the microstructures generated by the unit cells 
depicted in Figures j^a) and (b) are not the same. However, if the radii of the fibers and 
the spacing of the fibers relative to the width of the unit cell are chosen appropriately, then 
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Figure 5: The unit cell Y = ( — 1,1)^ consisting of a matrix embedded with fibers in different 
configurations. If the matrix and the fibers have tetragonal symmetry, then the associated 
macroscopic elasticity tensor has tetragonal symmetry, (a) The unit cell contains one fiber 
whose axis aligns with the y^-axis. (b) The unit cell contains four fibers oriented in the 
direction of the y^-axis. (c) The unit cell contains two fibers that are oriented in orthogonal 
directions. 


they do generate the same microstructure. As long as the material symmetry groups of 
the matrix and fibers contain and Re/^, the result still holds. Thus, if the fibers are 
transversely isotropic and the matrix is isotropic, then the macroscopic elasticity tensor 
has tetragonal symmetry. This type of periodic elastic structure was used by Ptashnyk 
and Seguin m in modeling plant cell walls. The macroscopic elasticity tensor also has 
tetragonal symmetry when the fibers are oriented in orthogonal directions, as depicted in 
Figure l^c). To see this, let 

{{yi,y2,y3) g T | (yi - 0.5)^ + y| < 0.2^}, 

{(yi, 2/2,2/3) G T I (?/i + 0 . 5)2 + 2/2 < q_22| 

be the domains occupied by the two fibers in the unit cell. It can be verified that the 
transformations 

z I—;■ Hz for all 2; G M^, H = —Rgj, “R-es 

and 

2; I-)- ei + Re(^z for all 2; G 

are symmetries of the elastic structure and, thus, the macroscopic elasticity has tetragonal 
symmetry. 
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Hexagonal symmetry: Let i? be a rhombus in whose short diagonal and sides have 
length 2. Consider the unit cell Y = ( — 1,1) x R. As in the example involving tetragonal 
symmetry, assume that Y consists of a matrix and hbers and that the hbers occupy the 
domain 


{(yi) 2 / 2 , 2 / 3 ) £ Y I dist((y 2 , 2 / 3 ),-R*) < 1, where Ri are the corners of i?} (23) 

and the rest of Y is occupied by the matrix. A cross-section of 1" at a constant yi value is 
shown in Figure [^a). The elastic structure resulting from tiling with this unit cell has 
hbers arranged in a hexagonal bundle structure; see Figure |^b). When the matrix and the 
hbers are isotropic, then the periodic elastic structure generated by this unit cell has the 
symmetry 

z I—)■ Zo + — Zo) for all z G M^, 

where Zo is any point on the axis of one of the hbers. It follows that the corresponding 
macroscopic elasticity tensor has hexagonal symmetry. 

TT / 3 

It is known that the form of an elasticity tensor possessing the material symmetry Re^ 
is the same as that with transverse isotropy; see Curtin [8]. Thus, the material symmetry 
group of the macroscopic elasticity tensor also contains for all 9. However, there is no 
mapping h of the form ^ whose gradient is R^^ when 9 is not a multiple of vr/S. Thus, 
this example shows that not every material symmetry of the macroscopic elasticity tensor 
is generated by a symmetry of the periodic elastic structure. 

In-plane fluidity: Suppose that the unit cell Y = (—1,1)^ consists of two linearly elastic 
materials, one occupying the bottom half (— 1 , 1 )^ x (— 1 , 0 ) ofT and the other the top half 
(— 1,1)^ X (0,1) of Y. Moreover, assume that both materials have in-plane fluidity in the 
plane orthogonal to 63 , meaning that their symmetry group is 

0 = {H unimodular | Hes = 63 and if 63 • a = 0, then 63 • Ha = 0}. 

This symmetry group is useful in the modeling of lipid bilayers; see, for example, Deseri, 
Piccioni, and Zurlo and Maleki, Seguin, and Fried |T2]. It is readily verihed that for all 
H G Q, the transformation 

z t-)- Hz for all 2 ; G 

is a symmetry of the periodic elastic structure and, hence. Theorem |3.2| implies that the 
macroscopic elasticity tensor also has in-plane fluidity. 

5 Conclusion 

Ill this paper we showed a connection between the symmetries of a periodic linearly elastic 
structure and the material symmetry of the macroscopic elasticity tensor obtained by using 
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Figure 6: (a) A cross-section with constant yi value of the unit cell Y = ( — 1,1) xi?, where R 
is a rhombus in whose short diagonal and sides have length 2. The grey regions, defined 
in (23), are occupied by fibers and the white region is occupied by a matrix, (b) A depiction 
of the periodic fiber structure obtained by tiling with the unit cell Y = (—1,1) x i?. The 
fibers, which are shown, are arranged in a hexagonal bundle structure. The surrounding 
matrix is not shown. Here ei is in the horizontal direction. 


homogenization theory. Our result generalizes previous work by considering a larger class 
of symmetries on the microscopic scale. Since the material symmetry group is a subgroup 
of the unimodular group, we conjecture that symmetries of the form (|^ include all pos¬ 
sible symmetries of a periodic elastic structure that lead to material symmetries of the 
macroscopic elasticity tensor. However, as the example in Section involving hexagonal 
symmetry shows, not all material symmetries of the macroscopic elastic tensor are generated 
by symmetries of the periodic elastic structure. 


The proof of Theorem 3.2 relies on the formula (12) for the macroscopic elasticity 


tensor. Since homogenization of problems posed in domains with locally-periodic or random 
microstructures also results in explicit formulas for the macroscopic elasticity tensor, it is 
possible that with suitable modifications the arguments used here could yield a result similar 
to Theorem 3.2 for these kinds of microstructures. It would be interesting to investigate 


whether there is an analog of Theorem 3.2 in the case of nonlinear elasticity. 


6 Appendix 

In this Appendix we present the details of the derivation of the macroscopic elasticity tensor 


given in (12). 


For the derivation, it is assumed that 

1. there is a strictly positive M such that |C(y)| < M for almost every y €Y, 
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2. the function y i—>■ C(y) is measurable, 

3. there is a strictly positive a such that a|Ep < E • CE for all E G Sym, 

4. A • CB = A • C[^(B + B^)] = [^(A + A^)] • CB for all linear mappings A and B, 
and 


5. D • CE = E • CD for all D, E G Sym. 

Let n be a reference configuration of the elastic material with elasticity tensor C^ defined 


in (111, and assume that is open and bounded with Lipschitz boundary so that locally the 


boundary can be represented as the graph of a Lipschitz function. Let Li and r 2 be disjoint 
subsets of with dO. = Li U r 2 , and assume a zero displacement boundary condition on 
Li and the traction on r 2 is given by t. The resulting mixed boundary-value problem in 
elastostatics with body force b is 


div(C=Vu^) -h b = 0 


in n. 


(24) 


= 0 on Ti, 

(C^Vu^)n = t on r2, 

where n is the exterior unit-normal to fl. Notice that by the assumed properties of C, only 


the symmetric part of Vu^ is relevant in (24). By the Lax-Milgram theorem HU a unique 


solution of (24) exists in the space 

= {u G V) I u = 0 on Li} 


(25) 


provided that b G L^(n,V) and t G L^(r 2 ,V). Moreover, the solutions of (24) are 
bounded in the Lf^-norm independent of e and, hence, there is a u G iL^(n, V) such that 
up to a subsequence 

converges weakly to u in V); (26) 

see HU. The goal of homogenization is to find an equation that characterizes u. 


To derive the macroscopic equations associated with (24), the notion of two-scale con¬ 


vergence first introduced by Nguetseng m and further developed by Allaire [2] is applied 
using the following function spaces. 


• Lpg^(y, V) is the set of all functions 'i/> from T to V that are T-periodic and Jy |'i/’(y)P dy 
is finite. 

• Hp^j.(Y,V) is the completion with respect to the iL^-norm of the space of smooth 
functions from T to V that are T-periodic. 

• Wper(E, V) consists of equivalence classes of V) where two functions are equiv¬ 

alent if they differ by a constant vector. 
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• L^(r2, Wper(^, V)) consists of all functions xp from Q, to Wper(^, V) such that 

In 

is finite. 

• C(n, Lpgj,(y, V)) consists of functions xp : Qx S ^ V such that for all x G 17 the func¬ 
tion y I—)■ xp{x,y) is in Lpgj.(y, V) and the function x i—)■ xp{x,-) from 17 to Lpgj,(y, V) 
is continuous. 


Two-scale convergence, in the context presented here, is defined as follows. 

Definition 6.1. A sequence of functions G L^(17, V) two-scale converges to G L^(17 x 

y,V) if for allxPeC{n,Ll^,{Y,V)), 

lim / (x) ■ xp(x, - — --\-q)dx= / f- 'v^(x,y) ■ xp(x,y)dydx, (27) 

e J^Jy 


where the symbol f denotes the average integral. 


By using the change of variables x x — q, one can see that this definition is equivalent to 
the standard definition mm involving a vector space rather than a Euclidean point space. 
Associated with two-scale convergence is the following compactness result 1211131 - 

Theorem 6.2. If G 77^(17, V) is a sequence of functions that weakly converges to G 
77^(17, V), then there zs a G L^(17, Wper(y, V)) such that up to a subsequence 


Vv*" converges two-scale to Vv^ -|- Vj^v^. 


It follows from Theorem 
that up to a subsequence 


6.2 


that (26) implies that there is a w G L^(17, Wper(y, V)) such 


Vu*" converges two-scale to Vu -|- V^yW. (28) 

Define v^(x) = v°(x) -|- S(p{x)v^{^^ q) for x G 17, where v° G 77 n C'°°(17,V), 
(p G and G 77pg^(y, V). Notice that 

v*" converges weakly to in 77^(17, V). (29) 


Multiply (24)i by v^, integrate over 17, and then integrate by parts to obtain 


— / Vu^(x) • C'^(x)[Vv°(x) -|- ev^(^^ — -q) iS) V(p{x)-\-(p{x)Vy 

Jn e 


1 .X - q 


+ q)]dx 


+ / t(x) • v'^(x) dx-|- / b(x) • v^(x) dx = 0. (30) 

Jr 2 Jn 
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Notice that the mapping 


(x, y) ^ C{y) [Vv°(x) + ^{x)VyV^{y)] 

is a suitable test function in the definition of two-scale convergence and that 

e / Vu^(x) • C^(x) ^-h g) <8) V(^(x)]dx—;> 0 as e —)• 0. (31) 

Jn ^ 

Taking the limit as e —)• 0 in (30) and using ( |28| ), (29), and (31) together with the symmetry 
of C yields 

- I f C{y)[Vu{x)+ Vyw{x,y)] ■ [Vv^{x) + 4>{x)VyV^{y)]dydx 
JnJY 

+ / t(x) • v°(x) dx-|- / b(x) • v°(x) dx = 0. (32) 

Jt 2 Jn 

Taking = 0 in (32) and using the arbitrariness of (/> G (^“(n) yields 

C{y) [Vu(x) -h Vyw(x, y)] ■ VyV^{y) dy = 0 

for all G ddpgj,(T, V) and almost every x G fl. (33) 

It can be shown that given Vu, there is exactly one w G L^(n, Wper(T, V)) that satisfies 
( |33[ ); see, for example, na. Also, w only depends on x through the symmetric part of 
Vu(x) and, moreover, w depends linearly on the symmetric part of Vu(x). Thus, (33) can 
be reformulated as follows: for any E G Sym, there is a unique G Wper(k^) V) such that 

^ C(y) [E + V 3 ,wE(y)] • V,v(y) dy = 0 for all v G V), (34) 

which is called the unit cell problem. 

Next, take </> = 0 in (32) to obtain 

- f C(y) [Vu(x)Vyw(x,y)] •Vv°(x)dydx 
JnJY 

-|- / t(x) • v°(x) dx -I- / b(x) • v°(x) dx = 0. (35) 

Jt 2 Jn 

If we define 

C°E '■= j~ ^(y) [E + Vj,w^(?/)] dy for all E G Sym, (36) 

which when written in components relative to an orthonormal basis has the form 


P n 

= f (Cijkiiy) + ^imiy)dy^^7®^'iy))dy^ 

p,q=l 


(37) 
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then (35) can be rewritten as 


— f C°Vu(x) • Vv°(x) dx + f t(x) • v°(x) dx + f h{x) ■ v^{x) dx = 0. 

J Q '^r2 


( 38 ) 


Since (38) holds for all G n C'°°(n, V), this equation is the weak formulation of 

div(C‘^Vu) + b = 0 in n, 

u = 0 on Fi, 

(C°Vu)n = t on r 2 , 


(39) 


which are the macroscopic equations associated with the microscopic system (24). This 
justihed calling the macroscopic elasticity tensor. A unique solution of (39) exists in T-L 
by the Lax-Milgram theorem na. Numerically, it is more efficient to solve the system ( |39[ ) 
than (24) since the elasticity tensor in (39) is constant, while in (24) the elasticity tensor 
can have rapid oscillations. 

An analogous homogenization argument can be carried out for elasticity problems in¬ 
volving different boundary conditions or including an inertial term. However, the resulting 
formula (12) for the macroscopic elasticity tensor, which is of primary interest here, would 
remain unchanged. 
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